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PART I 


SUMMARY OF WORK ACCOMPLISHED DURING THE PERIOD 
SEPTEMBER 1, 1971 - APRIL 30, 1972 

During the period of September 1, 1971 through April 30, 1972, 
several improvements were made on the Compensator Improvement Program 
(CIP). 1 In the first interim period (September 1, 1971 - December 31, 
1971) CIP was modified so that multi-channel systems and factored compen- 
sators could be handled. The factored compensators make allowances for 
the control of compensator pole-zero locations, and the factored compen- 
sators can be used to increase the convergence rate of the program. These 
two advantages of the factored compensators were demonstrated experi- 
mentally by using the CIP program . to design a single channel compensator 
for a Saturn launch vehicle that had been previously designed using a 
single rational function form. 2 

In the second interim period (January 1, 1972 - April 30, 1972) 
finalization of the CIP was accomplished. In this time interval several 
additional options were included in the program. These options are the 
choice of allowing nonminimum phase compensators, the choice of perturbing 
some frequency points with respect to more than one point, and the choice 
of the program continuance mode, i. e., the total improved frequency 
response mode or the sum improved frequency response mode. The practical 
usefulness of CIP was demonstrated by designing compensators for the 
Saturn V/Sl-C Dry Work. Shop and the Saturn V/Sl-C Sky Lab. From these 
examples it is pointed out that using CIP considerable improvement of 
compensators can be made, and the amount of effort to design compensators 
is reduced greatly. 3 
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As. a conclusion to CIP summaries of the subprograms were documented 
in May, 1972. This information was considered necessary in case some of 
the subprograms are to be modified or replaced by equivalent programs. 
These summaries present and explain the basic programming theory of the 
specialized subprograms (programs not normally found in a computing 
center system library). Also, complete descriptions of each I/O variable 
of each subprogram are given. With the programming theory and I/O vari- 
able information the modification or reproduction of any of these programs 
should be minimized. 4 
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PART II 

MAGNITUDE ERROR BOUNDS FOR SAMPLED-DATA FREQUENCY RESPONSE 
OBTAINED FROM THE TRUNCATION OF AN INFINITE SERIES 

Introduction 

In order to properly design a control system it is first necessary to 
be able to analyze the performance. In general there are three basic classi- 
cal approaches for analyzing control systems; these approaches are the 
frequency response method, the time domain method, and the root locus method. 
These methods all transmit the same information. However, deciphering some 
information which is obvious while using one technique can be very difficult 
while using either of the other techniques. Thus, each method has certain 
advantages and disadvantages which usually depend upon the type of system 
performance to be measured. 

Probably the most widely used technique of analyzing control system 
performance is by the frequency response method. The reasons that this 
technique is so popular are the simplicity of obtaining frequency response 
information with the aid of the digital computer and the availability of 
experimental frequency response information. 

The frequency response method has been of great utility in the analysis 
of sampled-data control systems that are unity feedback and error sampled. 
However, one major problem that is associated with the application of the 
frequency response method to sampled-data control system analysis is the 
difficulty of obtaining the sampled frequency response of the continuous 
system. A technique for accomplishing this is given by the following 
infinite series: 


G*(jco) = ^ l G(jm + jqu g ) + g(0)/2 


q — 0 


1 

T 


( 1 ) 
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where G ( j o>) is the sampled frequency response, G(jco) is the frequency 
response of the continuous system, T is the sampling period, u> s is the 
sampling frequency and o) is the radian frequency. 1 Practical use of (1) 
is made by truncating the series after a certain number positive and 
negative terms, e.g. , q f s on the interval of -k q k where k is a 
positive integer. The result of performing this produces the following 
finite series: 

1 k 

G*(jw) ® — l G( jw + jqw s ) + g(0)/2 (2) 

q=-k 

in which G*(juj) is the approximate sampled frequency response. It then 
follows that 

G* (jw) 2. G*(ju) . (3) 


Upper Bound Magnitude Error Derivation 


It is important to know how many positive and negative terms of the 
series must be taken so that G*(ja)) is a reasonable approximation to G*(jio) 

^S < (jj < ^ 

for - — _ . One technique for studying this is to determine the 

magnitude error between the two frequency responses. This error is defined 
as 


E(jto) 


G*(jco) - G*(j(o) 


(4) 


where E(joj) is the error between the actual sampled frequency response 
and the approximate sampled frequency response. Subtracting (2) from (1), 

~~ ju Cl) g CO g 

It is only necessary to consider the frequency range — — <_ uj , 

which defines t?he primary strip in the G* (s) -plane, since the 
frequency response G*(jio) is simply duplicated for frequencies in each 
secondary strip of the G*(s)-plane. 
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(4) becomes 
| E(joj) | 


00 —CO 

1 G(j to + jqu s ) +. I G (ju + jqoig) 

q=k+l q=-k-l 


(5a) 


or 


E(jw) 


1 

T 


l 

q=k+l 


[G (jto + jqw g ) + G(jto - jqu) g )] 


(5-b) 


In order to determine the magnitude error exactly it is necessary to 
evaluate the infinite series given by (5a) or (5b). This is impractical 
and actually defeats the purpose for truncating the original series. 
However, if an easy approach to calculating a reasonable upper bound on 
the magnitude error can be established, then a useful relation exists. 

Consider the continuous frequency response G(jco) and suppose that 
for all u) >_ k(jJ g and for all co <_ -ku) g that 

|G(jw)| £ M/|w n | (6) 

where M is an appropriately chosen constant and n is the difference in 
the number of poles and zeros of G(s). This assumption is certainly 
not unreasonable if ku) g is greater than all break frequencies of G(jto). 
Next assume that n >_ 1. 

Using the identity for complex numbers, z^ 9 z 2 , ... , zj which 
states that |zj + z 2 + ••• + z j | <_ | z ^ | + | z 2 | + ••• + | z j | (5) becomes 


* n 

It should be obvious that (6) assumes G(jo>) ^M/(jw) for all 

a) > ko) _ and for all u) < -kw . 

— o — * S 
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E(jca) 


1 

T 


l [G(ju> + jqu s ) + G(jio - jqa) )] 
q=k+l 


i ? l 


q=k+l 


G(ja) + jqw s ) + G(jw - jqoj s ) 


CO 

ill [ I G (jw + jqu s ) | + |G(ju) - jqo» s )|] 


q=k+l 


( 7 ) 


Using ( 6 ), the magnitude error is bounded by 


r-\l < M V r 1 , V 1 1 

J “ - T q= k+i L (<l“s + “ )n (< * w s ■ u > n J 




( 8 ) 


in which it has been assumed that qo) g > a) since. u)je: [- — , — ]. 

The result is that the magnitude error can be bounded by the sum of 
two infinite series. It is possible to bound either of these series from 
above and below by Riemann Zeta series minus certain constants. Since 
either of the series in (8) can be bound from below by a Riemann Zeta 
series then neither series converges for n=l because for n=l Riemann Zeta 
series diverge. Thus the above, bound is only applicable for system’s with 
at least two more poles than zeros, i. e., n >_ 2. 

Returning to (8) it is observed that both series are positive series. 
Thus the values of each series is less than the integral between the 


A Riemann Zeta series or p-series is an infinite series of the form 

OO T_ OO 1 

l/k. A series of the form ( qa)g + ^n where a> = pu) s and pe[p, 0 . 5 ] 

1 00 1 C ° ^ 1°°1 

is bound by — ^ ^ 1/ (qa) e + co) n > — ff ^ V • 

y w s n q=k+l q n - q =k+l 4 s “ q=k+l (q+D n 

The series on the extremes of this inequality are simply weighted Riemann 
Zeta Series with the first (k- 1 ) and k terms, respectively, subtracted off. 
The subtracted terms form the constants mentioned above. The series 

00 1 

Z. -7 3 — yn may be bound in a similar manner. 

q=k+l 
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limits of k and °° of the continuous curves fabricated by replacing 
each summing index by a continuous variable. 2 In mathematical notation 
this produces 


I E(joi). 


< E y r _i + J L 1 

~ T q=k+l L<^ S + w)n ^ U s " w)n J 


< *1 

- T 


- 

dx 

00 — 
dx 


j 

(xid + td) n 
k s * 

(xm g - u)n 

k J 


M 

T \ co s (-n+l) 


1 1 
(xoj + w) n- l + (xm - m) n_i 


(9a) 


or 


M 


| _ 2Tf(n-l) _ (ka) + w )n-l (kco - a)) n 


(9b) 


The expression given in (9) can be used to calculate an upper bound 

w s w s , . 

magnitude error for any ooe[- — , — ]• It is obvious that |E(jco)| is 

an even function with respect to co, i.e., | E ( jo>) | = |E(-ja))|. It is then 

only necessary to consider go’s on the range 0 co <_ — • 

' For co ^ s on this range it is easily deduced that the maximum value of 

the upper bound magnitude error occurs at to = co s /2. Thus the absolute 

maximum magnitude error that any point on this given frequency range can 

achieve is 

M 

E max 27ro) n “l(n-l) 
s 


(k + l/2) n_1 + (k - 1/2) 


n-!_ 


( 10 ) 


The relation (10) allows for the calculation of a bound on the maximum 
magnitude error that any point in the allowable frequency range. can 
possess. A quick estimate of the magnitude error of any frequency point. 
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where we[0, oj s / 2], can be obtained from (10), whereas (9) can be used to 
obtain a point by point maximum magnitude error for frequency points on 
the interval 0 co _< u> g /2. 

It is obvious from both (9) and (10) that the choice of M affects 
the maximum magnitude error upper bound. In fact the smaller M is, the 
smaller the upper bound error. Thus, the closer that M/|(o n j approximates 
the magnitude of G(ju)) for o> >_ ku) g and u) <_.-k(i> s the more accurate the 
upper bound magnitude error will be in predicting the actual error. One 
technique for calculating the value of the smallest M is to plot the 
| G ( j co) | for u> > ku3 g on log - log paper and then construct a straight line 
whose slope is [-(octave)x n/octave] which is greater than |G(jw)| for • 
u) > ku) g but is as near as possible to |G(jw) | . This straight line may be 
labeled as S(a)). From the log - log plot for some value of frequency, e,g. , 
co = co 2 > the corresponding value of S(u)) can be read, i.e., S(o)j). Then 
M = SO^) • 

Examples 

The previously derived upper bound magnitude error for bounding the 
magnitude error caused by truncating the infinite series^ (1), after k 
positive and negative terms is of practical use only if the bound can 
produce a reasonable estimate of the actual error. In order to illustrate 
that the upper bound calculation given in (9b) can do this , two examples 
of, its use are given in the followings 

Example 1 - Consider a continuous system described by the following 


G (s) 


1.0 

(s+l) 2 


trans f er function : 
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It is desired to use (9b) to approximate the sampled frequency response 
of this system when oj g = 25 rad/sec and k = 2. Table 1 is used to compare 
the approximate frequency response, to the exact frequency response and to 
compare the actual magnitude error to the upper bound magnitude error 
calculated from (9b) when M = 1.0.* 

Example 2 - In this case consider a continuous system described by 

(s+2) 

G s) = 256 ’° ( s + 4)(s + 8)(s + 16) 

For k = 5, u) s = 25 rad/ sec, and M = 256.0 Table 2 can be used to make 
comparisons between the approximate and the exact . frequency responses and 
the actual and the upper bound magnitude errors. 

Analysis of Results and Conclusion 

From the preceding examples it . is seen that the upper bound magnitude 
error equation, (9b) can yield reasonable estimates of the actual error. 
However, it should be pointed out that upper bound calculations depend 
very much upon M. In fact in the examples given the upper bound calcula- 
tions are as small as possible because M was chosen so that M/joj n | is the 
actual asymtote of jG(jto)| as to -*■ °°. If kw s had been small enough so that 
this would not have been possible then the upper bound calculation would 
have probably been much greater than the actual error, e.g., cases where 
ku> s is less than the last break frequency of G(jm) or the value of kw s is 
close to the break frequency of a pair of complex conjugate poles with 
6 < 0.707. 


Magnitude error is defined in (4). 
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There are two extensions or modifications which would greatly improve . 
upper bound magnitude calculations. First, the technique derived in this 
report is not applicable to systems with only one more pole than zero. 

Thus the upper bound magnitude derivation should be extended so that this 
case is included. Secondly, modifications of the derivation should be made 
so that a tighter bound can be calculated. This could be done by producing 
tables that could be used for the actual evaluation of (8) , This would 
eliminate the error caused by approximating the value of the infinite 
series of (8) by integrals of continuous variables. With these extensions 
the upper bound magnitude calculation would be enhanced greatly. 
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Table 2„ Frequency Response Information for Example 



